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Abstract 

We present a relation between the Mathieu equation and a particular elliptic 
curve. We find that the Floquet exponent of the Mathieu equation, for both 
q « 1 and q » 1, can be obtained from the integral of a differential one form 
along the two homology cycles of the elliptic curve. Certain higher order differ- 
ential operators are needed to generate the WKB expansion. We provide a fifth 
order proof. 
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1 Introduction 



Mathieu equation was first introduced by E. Mathieu when he studied vibrating elliptical 
membranes pp. Its canonical form is 

d 2 u 

— + (A - 2qcos2z)u = 0. (1) 
dz 2 

The related modified Mathieu equation is obtained by z — > iz: 

^4 - (A - 2gcosh2^)M = 0. (2) 
az 1 

The Mathieu equation is useful in various mathematics and physics problems. As an example, 
the separation of variables for the wave equation in the elliptical coordinates leads to the 
Mathieu equation. 

According to the Floquet theory, the solution of the Mathieu equation can be written in 
the form: 

u u (z) = e™ z f(z). (3) 

where f(z) is a function of period tt, and in general v is a constant independent of z. v 
is called the Floquet characteristic exponent, it is a function of the constants A and q. A 
classical result is that the Floquet exponent can be obtained through the Hill's determinant. 
Moreover, if v is an even integer, then the solution u(z) is a periodic function of period 7r; 
if v is an odd integer, then the solution u(z) is a periodic function of period 2n. In our 
discussion in this paper, u(z) is not required to be periodic. 

The Mathieu equation has been studied for a long time, for the collections of classical 
results see nice references [2J [3], HI [5] , and more recent studies in (6j [7] . 

Another object we study here is a particular elliptic curve. Geometrically the elliptic 
curve is topologically equivalent to a torus, it is a Riemann surface of genus g — 1. The 
relation between the Mathieu equation and the elliptic curve naturally aries in the integrable 
theory. The (modified) Mathieu equation is the Shrodinger equation of the two body Toda 
system, while the elliptic curve is just the spectral curve of the classical Toda system. See 
[8] for relevant backgrounds. As an illustration, let us start from the Mathieu operator 

C = d 2 z + A - 2gcos2z 

= rfJ + A-g^ + JL). (4) 

Substituting d z = x and q(e 2tz — e~ 2lz ) = y, where x,y are complex coordinates. Then we 
have 

£ = (x 2 + A) ± vV + 4g 2 . (5) 
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The relation 

y 2 = ( x 2 + A) 2 - 4q 2 (6) 

is nothing else but the elliptic curve we are interested in. 

The curve (jSJ) has two independent conjugate cycles a and j3, they are canonical basis of 
the homology class of the torus. According to the general theory of Riemann surfaces, there 
is a holomorphic differential one form on the torus: 

" = (7) 

y 

and we can construct two periods by integrating u along cycles a and (5. 

A=(pu, B = (hu. (8) 



Then r = 4, Imr > is the complex modulus of the elliptic curve. 
However, we are interested in a meromorphic one form, 

Co = — . (9) 

y 

It is related to w by u = — 2^ + J^(~)cfx, the total derivative term will not contribute to 
contour integrals. The reason for us to study u, rather than u, is that it is directly related 
to the Mathieu equation. As a first hint, let x 2 + A = 2gcos2z, then we have 



to 



y/X - 2qcos2zdz. (10) 



This is actually the leading WKB (Wentzel-Kramers-Brillouin) solution of the Mathieu equa- 
tion. In the next section, we will see that they have an even deeper connection. In physics 
literature, the elliptic curve (EJ) is called Seiberg-Witten curve, and u is the Seiberg-Witten 
differential [9]F . 

The elliptic curve ([6]) can be viewed as a double covering of the branched x-plane. There 



are four branch points at x — (iy/X + 2q, iy/X — 2q, —i\/X — 2q, —i\/X + 2q), and two branch 
cuts run between (iy/X + 2q, iy/X — 2q) and (—iy/X — 2q, —iy/X + 2q). The homology cycle a 
of the elliptic curve corresponds to the contour encircling singularities (iy/X + 2q, iy/X — 2q), 
and the homology cycle ft of the elliptic curve corresponds to the contour encircling singu- 
larities (iy/X — 2q, —iy/X — 2q). In the next two sections we will show that, for q « 1 the 



3 In some literatures, as in [9], the elliptic curve is presented in a cubic form y 2 = Yi^=i( x — e i)i while in 
some other literatures, as here, the elliptic curve is presented in a quartic form y 2 = J~[^ = i(x — e{). They are 
equivalent forms of the same curve. The cubic form is obtained from the quartic one by proper coordinates 
transformations and parameter redefinitions, the result is that one of the zeros is mapped to infinite. The 
cubic curve in [5] can be written as y 2 — {x 2 — fg)(x — -|), after identifying u = ^-A, A 2 = ^-q according to 
[TT] and rescaling x — >■ h 2 x,y — > H 3 y. 
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Floquet exponent v is given by integrals of differential one forms along the a cycle on the 
torus, for q » 1 the v is given by integrals of the same differential forms along the cycle. 

The relation between Mathieu equation and elliptic curve we present here is found in 
our study in[10j E], about a relation between gauge theories and quantization of integrable 
systems [T2]. It suggests us to develop a WKB formalism to solve the Mathieu equation, as 
we explain in the next section. In this paper we try to present the problem as a differential 
equation problem, for relevant physics background, see [12] and [TTJ H5] , and references 
therein. 



2 Floquet characteristic exponents from elliptic curve 

As the first step, we rewrite the Mathieu equation in a form convenient for WKB expansion. 
Suppose q » 1, we rewrite it as 

e 2 d 2 u 



2 dz 2 



+ (w - cos2z)u = 0, (11) 



where e 2 = ~,w = ^. Then e is a small expansion parameter. We expand u(z) as WKB 
series: 

/ \ i P p(z')dz i r (^^-+ Pl (z')+ep 2 (z')+---)dz' / 10 x 

Substituting the series expansion ( [121 into the equation (TTTT) , we can solve p(z) order by 
order. 

Of course, the requirement q » 1 is not always satisfied. One may wonder if the results 
we get can be applied to the case q « 1. As we will see later, by suitably adjusting A, 
we actually obtain two convergent series. One series is convergent for q » 1, -% << 1, 
surprisingly it is still valid for the region q << l,u » 1. Another series is convergent for 



The first few recursive relations for p m are: 



Po = y/2(w- cos22;), p 1 = U\np )' , 

P2 = --^[2(lnp )"-((lnp )') 2 ], Ps = 

8p 2 p 



(13) 



where the prime denotes J|. 

Then we extend the Mathieu equation and its periodic solution to the complex domain 
associated with the elliptic curve. Then p(z)dz is a differential one form associated to the 
elliptic curve. Actually, the leading order p (z)dz is proportional to the u we introduced 
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above. We are interested in the integrals of p(z)dz along the conjugate homology cycles 
a and (3 on the elliptic curve, or equivalently, along the contours encircling (— f , §) and 
(— |cos _1 w, TjCos^ 1 ^) on the z-plane. It is the monodromy of the Mathieu function along 
cycles a and (3 on the torus. 

The leading order integrals are related to the complete elliptic integrals of the first and 
second kind, the result is: 



<j)p (z)dz = ny/2(w + l)F{ 



TT w + V 



/ \ , ?7T/ , „,1 1 1 —W, , 

p (z)dz = ->-l)F(-,-,2;^-). (14) 



As pi,p3 are total derivatives, the contour integrals of them are all zero 

P2m+i{z)dz = 0, m = 0,l, (15) 



and 



a, /3 



^ 1 / sin 2 2z - Awcos2z + A ^ 

a / 2 Z ~ 8V5j a ,p {W ~ COS22) 5 / 2 

1 f cos22; 



-dz 



12V2 Ta,b {w - cos2z) 3 / 2 
= —(2wd 2 w + d w )({) y/2(w - cos2z)dz, (16) 

where d w — -4-. We have simplified the integral by discarding some total derivative terms, 
this method was first used in |13j . In a similar way we find 

1 28 8 5 f 

Vdz = ^(^ 2 di + + 3O j) Podz. (17) 

We can proceed the same technique to obtain the differential operators for higher order 
p m , by discarding total derivative terms and simplifying the expression as far as possible. 
We call these differential operators generating differential operators. Acting these differential 
operators on (ftpodz, we can get higher order contour integrals, they can be written as 
combinations of the hypergeometric functions by using the formula: 

4~F(a, b, c; z) = —F(a + 1, b + 1, c + 1; z). (18) 

dz c 

As a demonstration, the expression for j}p 2 dz can be found in [TT], and ^p^dz is even more 
lengthy. We can get series expansions near a suitable value of w from these hypergeometric 
functions. However, it is much simpler to get the series expansion of po first and then to act 
the generating differential operators on this series. 



4 



In principle, all higher order generating differential operators can be determined by WKB 
relations. However it turns out that the calculations become very involved and it is hard to 
determine whether the expressions can be simplified further by discarding total derivative 
terms. Based on some observation on Po,Pi,P2,P3,Pa, we make a conjecture for higher order 
differential operators. 

Claim 1: In general we have 

P2m+\dz = 0, 

J) P2 m dz = (c m , m w m di m + c m ^ l w m - l d 2 ™~ 1 ^ Po dz, (19) 

Ja,/3 Ja,/3 

where m — 0, 1, 2, • • • , and c m i , (i — 0, 1, • • • , m) are numerical coefficients. 

Now we will state the relation between the monodromy of the Mathieu function along a, (3 
and its Floquet exponent. The asymptotic expansions of hypergeometric function F(a, b, c; z) 
are quite different for z — 0,1, oo. For example, let us look at the asymptotic behavior of 
the leading order results <f aj3 p dz. At w = oo, we have 

v v ; v 2' 2' 'w + V [ 4 v 2w 64 v 2w 256 V 2w J ' 

X 111 id 1 11 

_(u,-l)F(-,-,2;-^) = -v^[(ln2«;-2 + 21n2) + -(1 - 21n2 - ln2 W )(— ) 2 

+ Y^(47 - 601n2 - 301n2 M )(i-) 4 + •••]. (20) 
While at w ~ 1, with a — w — 1, we have 

^T )F( J i,! * ) = i + ' 1 + 21 " 2 --'%' 3 - 41 "^ 21 "^ 

2 2 W + 1 7T 27T 047T 

3(2 -21n2 + lna)(7 3 

-|_ . . . 



5127T 

1. , „,1 1 1 - w, 1 1 9 3n 25 , .. 

-w-lF-,-,2: = -a a 2 + a 3 a 4 + • • • . 21 

2 V ; V 2' 2' ' 2 ; 2 32 512 16384 V ; 

It turns out that the asymptotic expansions of <§> a „ pdz which are only powers of w or a 
are related to the Floquet exponent of the Mathieu equation. 

Claim 2: The contour integral of p(z) along the a-cycle gives the Floquet 
exponent 

u = -(h p{z)dz, (22) 

fa 

for the case q » 1, 4j << I (or q « l,v » I), the hypergeometric functions 
should be expanded near A >> q » 1, i.e. w~oo. 
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Claim 3: The contour integral of p(z) along the (3 -cycle gives the Floquet 
exponent 

v = —<£p( z )dz, (23) 

2 

for the case q » 1,— << 1, the hypergeometric functions should be expanded 
near A ~ 2q, i.e. w ~ 1. 

In this way, we can get the function v = u(w, e) as series expansion of e and w. In order 
to obtain the eigenvalue A, we need to reverse the function v = u(w, e) to get w = w(u, e). 

3 5th order proof 

In order to prove the validity of our claims, we have to show that the asymptotic expansions 
of v given by the contour integrals are indeed the same as results known in literatures. This 
has been successfully done in [TT] for the first three orders e -1 po + ^Vi + e3 P4- m this section, 
we will show how to determine the generating differential operators of p^ and p$, following 
the Claim 1,2,3, which would be very involved for manual calculation. 
Let us start from a classical result of the asymptotic expansion for \ v : 

\ - 2 , 1 2 | 5z/ 2 + 7 4 

v - + 2{v*-l) q + 32(z/2 _ i)3 (z/ 2 _ 4)9 

9 ^ + 58 ^ + 29 *■ + ... (24) 



64(^ 2 - l) 5 (z/ 2 -4)(z/2-9) 

It often states in the literature that this asymptotic expansion is valid for q << 1 and v > 4. 
Actually, it is also valid in the parameter region q » 1 and << 1, this makes our WKB 
method applicable. Then we reverse the series (l24"j) to obtain the series for v as a function 
of A, q. This can be easily achieved with the help of computer programs, for example the 
Mathematica software. We can trust the inverse results up to the order q 6 . 
The inverse series gives 

v = v/a-^A- 3 / 2 -^- 5 / 2 
4 4 

-7/2 _ _|_ 35g 4 ^_ 9/2 



v 4 


+ 


15 «V 

64 ; 




+ 


273g 4 


1 4 


64 




+ 


4147g 4 


v 4 


64 




+ 


65637g' 


v 4 


64 



4 32 

105g 6 11/2 g 2 33g 4 1155g 6 

256 J 1 4 2 256 ' 

5005^ _ «2 8229g 4 42185g 6 17/2 

128 J 1 4 32 128 J 

4 72200V u -i9/2 ,<? 2 ^ 65569g 4 6294301^ 



(9(A" 23/2 ). (25) 
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We have cut off the A expansion at 0(X 23 / 2 ), and discarded all the q expansion terms 
beyond the scope of the accuracy of (1241) . Rewrite the inverse series in w and e: 



v = i[(2 W )^-^( 2w )-3/ 2 _^ (2w) -V2_l| (2w) -n/ 2] 
+ e [-^)- 5/2 - |(2^)- 9/2 " (2u0- 13/2 ] 

+ ^[-\{2w)-^ - | (2w)-^ - ^(2w)-^] 

+e 7 [- 1 -(2w)- U/2 - ^V)~ 15/2 " 7 -^(2w)-^]. (26) 

It is straightforward to expand the integral <§> (e _1 p + ep 2 + e 3 P4)dz at w = oo, and compare 
with (J26]). They indeed match [TUl E] • 

Interestingly, the expansion ( 1261) is precise enough to further determine the differential 
operator for e 5 p^. According to our first claim, we set 

d) p 6 dz = {c^ 3 w 3 d 6 w + c 3}2 w 2 d 5 w + c 3A wd A w + c 3fi d 3 w ) (h p dz. (27) 

Ja,f3 Ja,p 

Expanding (f) a p dz = a/2(u> + 1)F(— 1/2, 1/2, 1; 2/(w+ 1)) at w = oo as in ([2D]) , and acting 
the 6th order differential operator (1271) on the series, we get 



j> Pedz = - ^(315c 3 ,3 - 70c 3 ,2 + 20c 3 ,i - 8c 3 , )(2w;)- 5/2 
105 

- (1287c 3 , 3 - 198c 3 , 2 + 36c 3il - 8c 3 , )(2u;)- 9 / 2 



32 

- ^^(3315c 3i3 - 390c 3i2 + 52c 3>1 - 8c 3 , )(2 W )- 13 / 2 
225225 

- ^7^-( 6783c 3,3 - 646c 3 , 2 + 68c 3>1 - 8c 3 , )(2«;)- 17 / 2 

TOT /I 7^7^ 

- 131Q72 (12075c 3 , 3 - 966c 3 , 2 + 84c 3jl - 8c 3 , )(2^)- 21 / 2 

+ .... (28) 

In order to determine the four coefficients c 3j j, we have to match (128"]) with terms of order e 5 
in formula ( 126]) . A crucial point is that although there are only three nonzero terms in (126]) . 
limited by the accuracy of ( 12~4"]) to the q 6 order, the leading term w~ 5 ^ 2 is absent in ( 12 6p . this 
fact enables us to determine the four coefficients in (1281) . We finally arrive at 

1.124 3j6 158 2j5 153 j4 41 j3 , I 
^ = 2^ ( 945 ' + 105 w + ~35 w + 14 w) f ( ^ 

Terms of order w~ 21 / 2 and higher in ( 128]) are superfluous for the determination of c 3j j. After 
c 3 i are determined, these higher order terms can be subsequently determined, too. 
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The e 7 ps order contour integral has five unknown coefficients 

G) p s dz = (c 4A w 4 d 8 w + c^ 3 w 3 d 7 w + c^ 2 w 2 d 6 w + c 4A wd 5 w + c 4i0 cZ 4 ) (h p dz. (30) 

In order to determine the coefficients we need at least four terms in the e 7 order in (1261) . as 
the coefficient for e 7 {2w)~ 7 ^ 2 should vanish. Therefore, we need the g 8 order contribution for 
X u . Fortunately, it has been worked out in [7J: 

1469z/ 10 + 9144z/ 8 - 140354z/ 6 + 64228i/ 4 + 827565z/ 2 + 274748 „ 



I192(v 2 - 16)(u 2 - 9)(z/ 2 - 4) 3 (z/ 2 - l) 7 



It extends the e order terms in fl26l) to 



e 7 [- L -{2w)-^ 



-Q 



4147 ^ _ I!!007 (2lt , rl8 / 2 _ 1000684685 ^ 



64 



256 



16384 



This determines e 7 ^ ps^- 2 as 

1, 127 4 „ 13 



«,/8 



ndZ = 2^4725^ dw+ T75 W d " + 63^¥ W dw+ 945^¥ Wd - 



15229 4 
135 x 2 7 ™ 



(31) 



(32) 



Podz. 
(33) 



We have shown that the Claim 1 is correct up to the 5th order. Moreover, by using 
Claim 2 for q << 1, we have determined all the coefficients in the generating differential 
operators of p^ and p$- Then it is straightforward to expand §^pdz to the e 7 order, near 
w ~ 1, to obtain the Floquet index z/ = u(w,e) for g >> 1. After reverse the series 
v = u(w, e) to w = w(z/, e) and rewrite it in v, A, q, we get 

4^2 _ l 4z/ 3 - 3u 



Xu 



2 *"W*-. 2 3 ■ 26 ^ 
80z/ 4 - 136z/ 2 + 9 528z/ 5 - 1640z/ 3 + 405z/ 



2 12 g 2«gi 




2016z/ 6 - 10080z/ 4 + 5886z/2 - 243 




' 2 19 g 2 




33728z/ 7 - 249872z/ 5 + 276004z/ 3 - 41607z/ 




224gf 




2403072z/ 8 - 24881920z/ 6 + 45534368z/ 4 - 160£ 


^7536z/2 + 506979 


' 2 31 g 3 




44811520z/ 9 - 620967168z/ 7 + 1724770656z/ 5 - 

„„ 7 


1152647184z/ 3 + 130610637// 



2™qi 



(34) 



We keep only terms consistent with the accuracy limit e 7 ps- This is the classical result of 
the Mathieu equation for q » l(See formula (20.2.30) in [5J). This finishes the proof of our 
Claim 3. 
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4 Concluding remarks 



We show that the Mathieu equation is closely related to an elliptic curve, therefore there 
is a geometric structure for the Mathieu equation which is not captured by asymptotic 
analysis. The Floquet exponent of the Mathieu equation can be derived from integrals of 
certain differential forms along the homology cycles of the curve. These differential forms are 
determined by the WKB procedure. Integrals along each homology cycle give an asymptotic 
expression for the Floquet exponent expanded at a specific point, the inverse series gives 
the corresponding eigenvalue. Integrals along all homology cycles a and give the complete 
asymptotic expansions ff2T| and flMl) for the eigenvalue. 

The appearance of Riemann surfaces associated with differential equations is quite fa- 
miliar in the theory of integrable models, the Riemann surfaces are the spectral curves of 
the classical integrable system while the differential equations are the Shrodinger equation of 
the same system. The (modified) Mathieu equation we discuss here is simply the two body 
Toda system, the eigenvalue formulae ( 1241) and ( 13T1) are states with large quantum num- 
bers(expanded at A >> 1) and states with small quantum numbers (expanded at A ~ 2q), 
respectively|10[ ITT] . In fact, the Lax matrix L(z) of the classical two body periodic Toda 
system is[TT] 

rf \ - ( Pi e m-» + z -i e »-ffi \ 

L[Z) ~ \ + ze" p 2 J [6b) 

where P\,P2 are momentum and (71,(72 are coordinates of the particles, z is the complex 
spectral parameter. The spectral curve is obtained by the determinant 

det(/zI-L(z)) = 0. (36) 

It's of the form 

2t(ji) -(z + z- 1 ) = 0, (37) 

where 2t(/i) = ji 2 — (pi + p 2 )fi + f(qi — q2,PxP2) is a polynomial of degree two. By a change 
of variables, x = p,,y = z — £(//), the curve becomes 

y 2 = t 2 (x) - 1. (38) 

This is the same curve as (jSJ) modula a further coordinates change. Now let us substitute \x 
by the differential d z , then formally we can get a second order differential equation through 

det(W z - L(z)). (39) 

It's easy to see that after some proper coordinates changes this is the Mathieu differential 
operator. So in some sense, the Mathieu equation defines a "quantum torus". Actually, 
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Talalaev has constructed the "quantum spectral curve" for general integrable systems along 
this way [13]. 

It is possible that the relation we present here is just a particular case of a general pic- 
ture. As an example, the spectral curve for the two body elliptic Calogero-Moser integrable 
system is an elliptic curve closely related to ([6]) (elliptic curves always can be written in the 
Weierstrass form), and its quantization leads to the Lame equation. See a recent discussion 
in [16]. 

Another example is the direct generalization of the case we present here, the spectral 
curve for the iV-body An periodic Toda chain is a hyperelliptic curve of genus g = N — 1, 
with N > 3. The Gutzwiller's quantization scheme of periodic Toda chain introduces a rather 
involved Bethe-like quantization condition which involves both the Floquet exponents and 
the integrals of mot ion [TT]. Recently in [18] the quantization condition has been rewritten in a 
functional form that only involves the Floquet exponents of the associated Hill's determinant. 
We wonder if the idea presented here can be generalized to higher genus curves and provide 
a solution to the eigenvalue problem of periodic Toda chain. Note that for a iV-particle Toda 
system, there are N — l independent Floquet exponents Ui, i — 1, 2, • - • ,N — l(the condition 
^2f =1 v% = just reduces the center of mass motion, or is the traceless condition for the SU(N) 
group), and for the associated hyperelliptic curve there are 2(iV — 1) independent homology 
cycles aj and There are also a meromorphic one form and its WKB descendants on the 
curve. The differential forms involve exactly N—l coefficients Ik, with k = 2,3, ■ ■ ■ ,N, which 
are the integrals of motion of Toda chain. Among the integrals of motions I 2 is interpreted 
as energy while Ik for k > 3 have no physical interpretation. There is evidence that the 
Floquet exponents z/j are given by integrals of these differential forms along the homology 
cycles on the hyperelliptic curve [19]. This is enough to determine the functional relations 
between the Floquet exponents and the integrals of motion z/j = Vi{I 2 ,I^,- •• ,In)- There 
are 2(N — 1) asymptotic expansion points at Ik » 1 and at the dual points Ik ~ 1^ ■ The 
critical values lj® are determined by the Chebyshev polynomial [20] • If the N — l integral 
cycles are chosen as («{«}, /%}), satisfying {i} C {1,2, •••iV — 1}, {j} C {1,2, •••AT — 1} 
and {i} U {j} = {1, 2, • • • ,N — 1}, {i} n {j} = 0, then functions v {i) = v {i} {I 2 , h, ■ • • , I N ) 
have asymptotic expansions at Ik » 1, while vyy = V{jy(I 2 , 13, ■ ■ ■ ,In) have asymptotic 
expansions at Ik — 1^ << 1. By reversing the Floquet exponents Ui = Vi{I 2 , h, ■ • • , In) we 
obtain the eigenvalues Ik = hi^i, v 2 , ■ ■ ■ , fjv-i)- 
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